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Abstract 



(— I ' We study the mfluence of the chiral phase transition on the chiral magnetic effect. The chiral 

electric current density along the magnetic field, the electric charge difference between on each 
! side of the reaction plane, and the azimuthal charged-particle correlations as functions of the 

' temperature during the QCD phase transitions are calculated. It is found that with the decrease of 

the temperature, the chiral electric current density, the electric charge difference, and the azimuthal 
, charged-particle correlations all get a sudden suppression at the critical temperature of the chiral 

phase transition, because the large quark constituent mass in the chiral symmetry broken phase 
^ quite suppresses the axial anomaly and the chiral magnetic effect. We suggest that the azimuthal 

■ charged-particle correlations (including the correlators divided by the total multiplicity of produced 

charged particles which are used in current experiments and another kind of correlators not divided 
by the total multiplicity) can be employed to identify the occurrence of the QCD phase transitions 
in RHIC energy scan experiments. 
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I. INTRODUCTION 



The phase transitions of quantum chromodynamics (QCD), for example the phase tran- 
sition of the chiral symmetry restoration and the deconfinement phase transition, have at- 
tracted lots of attentions in recent years. It is expected that these phase transitions occur 
and the deconfined quark gluon plasma (QGP) is formed in ultrarelativistic heavy-ion col- 



lisions 



[liJs 



lider 



5| (for example the current experiments at the Relativistic Heavy Ion Co 
(RHIC) and the Large Hadron Collider (LHC)) and in the interior of neutron stars 
Furthermore, studying the QCD phase transitions is also an elementary problem in strong 
interaction physics. 

Recently, The STAR Collaboration at RHIC report their measurements of azimuthal 
charged-particle correlations near center-of-mass rapidity in Au + Au and Cu + Cu collisions 
at a/stvtv = 200 GeV. They find a significant signal consistent with the charge separation 
of quarks along the system's orbital angular momentum axis [o], [lo|. The observed charge 
separation indicates that parity-odd domains, where the parity (V) symmetry is locally 



violated, might be created during the relativistic heavy-ion collisions jUl, Il2 |. 

QCD is an SU(3) Yang-Mills gauge theory coupled with quarks. The gauge field can have 
nontrivial configurations which can be characterized by a topological invariant, the winding 
number [l3|. The winding number is an integer and reads Qw = -:§2^ J d'^xG'^^^G'^^ . 
Here g is the QCD coupling constant. The gluon field tensor and its dual are GJ^^^ and 
^'a' ~ \^^^''^G^pa- The nontrivial gauge field configurations with non-zero winding number 
Qw can result in non-conservation of the axial currents due to the axial anomaly i.e., 

d^jl = 2j2mf{^ft^,^Pf)A - ^G'^l.^r, (1) 
/ 

where ipf is a quark field, is the current mass of the quark, and Nj is the number of 
quark flavors, = J2f{i'f1fj.l5'^f)A denotes the axial current density in the background of 
a gauge field configuration A^. Integrating the two sides of Eq. ([T]) over three dimension 
space and one dimension time and assuming the number of right-handed and left-handed 
fermions is equal initially at t = — oo, we obtain 

{Nr - NL)t=oo = -2NfQ^. (2) 

From Eq. ([2]) one can clearly find that through the interactions between quarks and non- 
trivial gluon configuration with non-zero Q^j, the right-handed quarks are converted into 



left-handed quarks, vice versa depending on the sign of the winding number. Then there will 
be an asymmetry between the number of right- and left-handed quarks. When a magnetic 
field is added, an electric current is induced along the magnetic field and positive charges 



15|-|20|. 



are separated from negative charge, which is called the "chiral magnetic effect' 

Now that the chiral magnetic effect can be observed through the measurements of az- 
imuthal charged-particle correlations in the relativistic heavy-ion collisions, a natural ques- 
tion arises, i.e. whether can we detect the properties of the QCD phase transitions, especially 
the chiral phase transition through the observations of the chiral magnetic effect? To answer 
this question, we have to study how the chiral magnetic effect or the charge separation effect 
is influenced by the chiral phase transition. This is our central subject in this work. 

This work is the extension of our former work f2l| and presents many details. The paper 
is organized as follows. In Sec. II we simply introduce the thermodynamics of the 2+1 
flavor Polyakov-Nambu-Jona-Lasinio (PNJL) model. In Sec. Ill we will calculate the chiral 
electric current density along the direction of the magnetic field and study its dependence on 
the temperature during the QCD phase transitions. In Sec. IV we will calculate the electric 
difference between on each side of the reaction plane and the azimuthal charged-particle 
correlations in heavy ion collisions. In Sec. V we present our summary and conclusions. 



II. THERMODYNAMICS OF 2+1 FLAVOR PNJL MODEL 



In this work, we will study the chiral magnetic effect and the QCD phase transitions 
in the 2+1 flavor Polyakov-Nambu-Jona-Lasinio model. The validity of the PNJL model 
has been confirmed in a series of works by confronting the PNJL results with the lattice 



QCD data 



22h26| . The PNJL model not only has the chiral symmetry and the dynamical 



breaking mechanism of this symmetry, which are same as the conventional Nambu-Jona- 
Lasinio model, but also include the effect of color confinement through the Polyakov loop. 
Therefore, the PNJL model is very appropriate to describe the QCD phase transitions at 
finite temperature and/or density. 

In this work, we employ the 2+1 flavor Polyakov-loop improved NJL model which has 
been discussed in details in our previous work 25|], and the Lagrangian density for the 2+1 



3 



flavor PNJL model is given as 

8 
a=0 

+ {^il5ra^'f] - ir[det/(V^(l +75)^) 

+det;(^(l-75)^)] -W($,<l>%r), (3) 

where ip = {ipu, ipd, i^s)'^ is the three-flavor quark field, 

D'' = d''- iAf" with A" = , A° = c/^°y = -iA^. (4) 

Xa are the Gell-Mann matrices in color space and the gauge coupling g is combined with 
the SU(3) gauge field A'^{x) to define A^{x) for convenience, mo = diag(mQ, mg, mp) is the 
three-flavor current quark mass matrix. Throughout this work, we take rriQ = ttlq = rriQ, 
while keep mg being larger than mp, which breaks the SU{3)f symmetry. In the above 
PNJL Lagrangian, U ($,$*, T) is the Polyakov-loop effective potential, which is expressed 
in terms of the traced Polyakov-loop $ = {TtcL)/Nc and its conjugate $* = {TtcL'^)/Nc 
with the Polyakov-loop L being a matrix in color space given explicitly by 

If" ^1 

L{x) =V exp i I dr A4 (x, r) = exp [i/3A4] , (5) 

I Jo J 
with P = 1/T being the inverse of temperature and A4 = lA^. 

In our work, we use the Polyakov-loop effective potential which is a polynomial in $ and 
$* Q , given by 

+ (6) 

with 

&2(T) = ao + ai(|)+a.(^)%a3(|)'. (7) 

Parameters in the effective potential are fitted to reproduce the thermodynamical behavior 
of the pure-gauge QCD obtained from the lattice simulations. Their values are ao = 6.75, 
Oi = —1.95, 02 = 2.625, 03 = —7.44, 63 = 0.75 and 64 = 7.5. The parameter Tq is the critical 
temperature for the deconfinement phase transition to take place in the pure-gauge QCD and 
To is chosen to be 270 MeV according to the lattice calculations. Furthermore, we also need 
to determine the five parameters in the quark sector of the model, which are itLq = 5.5 MeV, 
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= 140.7 MeV, GK^ = 1.835, KA^ = 12.36 and A = 602.3 MeV. They are fixed by fitting 
= 135.0 MeV, tuk = 497.7 MeV, m^/ = 957.8 MeV and U = 92.4 MeV 
In the parity-odd domains which are created during relativistic heavy-ion colhsions, the 
number of left- and right-hand quarks is different because of the axial anomaly. In this work 
we introduce the chiral chemical potential /is to study the left-right asymmetry following the 
method of Ref. 15|, where the chiral chemical potential fi^ is related with the effective theta 



angle of the 6*- vacuum through /is = 3^6 /2Nf and Nf is the number of flavor. Consequently, 
we should add the following term 

#57075^ (8) 

to the Lagrangian density in Eq. ([3]), where /is = diag(/i5, /Xg, /ig). Next, we consider the case 
that a homogenous magnetic fleld B is along the direction of the orbital angular momentum 
of the system produced in a non-central heavy-ion collision. In the following we denote 
this direction with 2;-direction and particle momentum in this direction with ^3. We derive 
the corresponding thermodynamics for a system with only one kind of a fermion in detail in 
Appendix lAl In the Appendixes, we have emphasized that the approach used in Appendix lAl 
is appropriate to describe the chiral magnetic effect, while the modified Lagrangian approach 
in Appendix [B] is inappropriate. The results in Appendix |X] can be easily generalized to 
the 2+1 fiavor PNJL model and in the mean field approximation, the thermodynamical 
potential density (fi = —{T/V) InZ) for the 2+1 fiavor quark system under a homogeneous 
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background magnetic field B and with left-right asymmetry is given by 

f=u,d,s n=0 s=±l ^ 



+1 In { 1 + 3$* exp [ - [E, - - ^/.{) /T] 



-3<l>exp[-2(E^-/i^-^/iO/T] 



+ exp [ - 3 [Es -fif- '-^fil) /T] } + I In { 1 
+3$ exp [ - [Ef + fif- '^d)/T] 



+3$* exp [ - 2 {Ej + /i^ - /T] 

Ef 

+ exp[-3(E, + /.,-^/.9/T]}) 

+W($,$*,T), (9) 



where we have 



and 



ef\ = J2n\qf\eB+pl (10) 



Ef = ^2n\qf\eB + pl + M]. (11) 

Here = u, d, s) is the electric charge in unit of elementary charge e for the quark of flavor 
i. The constituent mass Mj is 

M, = mi - + 2K(Pj 0^, (12) 

and (pi is the chiral condensate We also include the quark chemical potential fii in 

Eq. (Q. From our calculations in Appendix [Aj one can find that the momenta of charged 
particles in the longitudinal direction, i.e., the z-direction, are not influenced by the back- 
ground magnetic field and in the expression of the thermodynamical potential density 
in Eq. is continuous; while the momenta in the transverse plane are discretized due to 
the external magnetic field. We should emphasize that at the lowest order of the transverse 
quantum number n = 0, i.e., the lowest order Landau level, the quark spin only has one 
value in the z-direction, which means that charged particles in the lowest order Landau level 



6 



are polarized by the external magnetic field; however particles in higher levels, i.e., n > 0, 
are not polarized. Therefore, the charge separation effect only comes from quarks in the 
lowest order Landau level. 



III. CHIRAL ELECTRIC CURRENT ALONG THE DIRECTION OF THE MAG- 
NETIC FIELD 



The chiral electric current along the longitudinal direction, i.e. the direction of the 
magnetic field, is an observable which describes the magnitude of the charge separation 
effect. Here, we use the approach in Appendix |A] to calculate the chiral electric current 
density js, whose expression is 



J3 



f=u,d. 

r>00 



qje'B 



Jo Ef hf Jo hf 



P3 „/ 



Pi II 



dps^f{Ef + f^f + ^f,i) 
Ef 



(13) 



Ef JO 

where qe = diag(g„e, g^e, g^e) is the electric charge matrix for three-flavor quarks; V is the 
volume of the system and Ef is given by Eq. ( ITTi) with n = 0. We have 



and 



fix) 



(14) 



(15) 



1 + 3$e-^/^ + 3$*e-2^/^ + e'^^/^ ' 
In order to study whether the chiral electric current is affected by the temperature, quark 
chemical potential, quark constituent mass and so on, we just pause here, and turn to the 
simpler system composed of only one type of fermion with positive charge e and mass m. 
The chiral electric current corresponding to this system is given by 

e^B \ r , P3 1 r . 1 



J3 



,0 '^^'Ee(^-^-t'^-)/^ + l 



dp: 



dp: 



E g(£;-M+^M5)/T _^ I 

P3 1 



(16) 
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where 

E = + m2. (17) 

It is interesting to consider the case that fermions are massless, i.e., m = in Eq. ( fTTl) . then 
it can be easily obtained that 



J3 



2n' 



(18) 



This is the result obtained in the modified Lagrangian approach as Eq. (]B23|) shows. How- 
ever, the chiral electric current calculated in our approach is in essence different from that 
obtained in the modified Lagrangian approach. From our calculations above, one can find 
that the chiral electric current comes from the finite temperature part of the thermody- 
namics, i.e., from fermions and anti- fermions, whereas the electric current in the modified 
Lagrangian approach comes from the Dirac Sea, not from fermions and anti-fermions (for 
more details see Appendix [B]). Furthermore, although obtained in our approach only de- 
pends on the magnetic field strength B and the chiral chemical potential /is in the massless 
case, it is indeed dependent on the temperature and the chemical potential fi when the mass 
of fermions is nonvanishing. On the contrary, in the modified approach only depends 
on B and fi^, regardless of whether the mass of fermions is vanishing, which is due to the 
unphysical ultraviolet momentum in the Dirac Sea. 
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FIG. 1: (color online). Left panel: chiral electric current density J3 in Eq. (I16p as function of 
the temperature with several values of the fermion mass m, with eB = iC^MeV^ /i = 0, and 
/is = 250 MeV. Right panel: js as function of the temperature with several values of the chemical 



potential /i and with eB = lO'^MeV^ /is = 250 MeV, and m = 200 MeV. 
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In order to verify our argument, we calculate the chiral electric current density in Eq. ( !T6|) 
for the one fermion system numerically, and the results are plotted in Fig. [1] From the figure 
one can easily find that, although the chiral electric current density is independent of the 
temperature when the fermion is massless as the black solid line in the left panel shows, it is 
indeed dependent of the temperature and the chemical potential in the case that the mass 
of fermion is nonvanishing. The effect of the mass of the electric current carrier, i.e., the 
fermion, is to decrease js, since the velocity of the fermion becomes smaller when the mass 
of the fermion is increased with a fixed energy. On the contrary, the chiral electric current 
density increases with the temperature and the chemical potential. 
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FIG. 2: (color online). Constituent masses of u, d quarks and s quarks, the Polyakov-loop <I> and 
its conjugate <I>* as functions of the temperature with = (i = u,d,s) and ^5 = 250 MeV 
(^5 = ^5 = ^5 = in the PNJL model with parameters given in the SecUll 



Next, we turn our attentions to the chiral electric current produced in the three- flavor 
quark system under an external magnetic field, whose expression is given in Eq. (fT3|) . From 
our above experience that the chiral electric current density would be affected by the particle 
mass and temperature, we first investigate the dependence of the constituent masses of 
three-flavor quarks on the temperature during the QCD phase transitions in the PNJL 
model. Minimizing the thermodynamical potential in Eq. with respect to three-flavor 
quark condensates, the Polyakov-loop $ and its conjugate $*, we obtain a set of equations 
of motion. We neglect the influence of the magnetic field on these equations of motion in our 



numerical calculations, since the mae: netic field {eB = 10^ ~ lO^MeV^ in the non-central 



heavy- ion collisions [1^) has little impact on these equations of motion. The calculated 
results are presented in Fig. [2l Here we take the chiral chemical potential fi^ = 250 MeV for 
example. In the figure we can find that the constituent masses of u, d quarks and s quarks 
decrease with the increase of the temperature, and a first order chiral phase transition takes 
place at the critical temperature Tc = 185 MeV. The chiral symmetry is restored above 
this critical temperature and the constituent masses of u, d quarks decrease to their small 
current masses. Since the s quark has relatively larger current quark mass, its constituent 
mass is still relatively large when the temperature is larger than T^, but it is also decreased 
quickly with the increase of the temperature. In our calculations we also note that when the 
chiral chemical potential /is is decreased, the value of the critical temperature Tq becomes 
larger and the first order chiral phase transition gradually evolves to a continuous crossover. 
Therefore, the influence of the chiral chemical potential /is on the chiral phase transition is 
similar with that of the quark chemical potential /i. In Fig. |2]we also plot the Polyakov-loop 
$ and its conjugate $* versus temperature. One can find that $ and $* increase from to 
1 with the increase of the temperature, implying that the Z(3) symmetry of the gluon field 



is broken and the deconfinement phase transition takes place |22l]. We should emphasize 
that since the term related with the chiral chemical potential in Eq. ([8]) does not broken the 
charge conjugation symmetry, we have $ = $* even /is is nonvanishing, which is different 
from the quark chemical potential. 

In Fig. [3] we show the chiral electric current density of the three-flavor quark system as 
function of the temperature in the PNJL model. Here we take the external magnetic field 
eB = lO^MeV^ for example. First of all, we consider the high temperature limit. In this 
limit the masses of quarks can be neglected and $ = $* = 1. Then the chiral electric current 

in Eq. f|T3|) can be easily obtained as 



f=u,d,s 



We also plot these high temperature limit values of the chiral electric current density in 
Fig. [31 i.e. the horizontal lines from bottom to top corresponding to /is = 150, 200, and 
250 MeVs, respectively. One can clearly find that when the temperature is high, the system 
is in the chiral symmetric and deconfined phase, and the chiral electric current density 
approaches its limit value, i.e. Eq. ( IT9|l . When the temperature is lowered, especially 
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FIG. 3: (color online). Chiral electric current density along the direction of an external magnetic 
field of the three-flavor quark system js as function of the temperature in the PNJL model with 
eB = lO^MeV^, = [i = u,d,s), and several values of the chiral chemical potential /U5 (//s = 
/i5 = /i5 = The three horizontal lines denote the values of J3 in the high temperature massless 
limit, corresponding to /^s = 150, 200, and 250 MeVs from bottom to top, respectively. 

when the temperature is below the critical temperature Tq, the chiral symmetry is broken 
and quarks get large constituent masses, then the chiral electric current density along the 
direction of the external magnetic field is quite suppressed and quickly approaches zero with 
the decrease of the temperature. This behavior is independent of the value of the chiral 
chemical potential as Fig. |3] clearly shows. In Fig. [3] we just take magnetic field strength 
eB = 10^ MeV^ for example, and the chiral electric current density is linearly proportional 
to the magnetic field strength, since the magnetic field strength B does not enter into 
the integrations in the expression of js, i.e., Eq. ( !T3l) . Therefore, when the magnetic field 
strength takes other values, we still have the fact that when the temperature is below the 
chiral critical temperature Tq, large constituent masses of quarks suppress the chiral electric 
current drastically. 
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IV. AZIMUTHAL CHARGED-PARTICLE CORRELATIONS IN HEAVY-ION 
COLLISIONS 



In this section we will try to relate our calculations with experimental observations and 
investigate how the QCD phase transitions influence on the signals of the chiral magnetic 
effect. In the experiments of heavy ion collisions, the azimuthal charged-particle correlations, 
i.e., (cos(0a + 0/3 — are used to detect the P-violating effect j9,[lo|, 28l. Here 0q, and 

0/3 are the azimuthal angles of the produced particles, and a, /3 represent electric charge + 
or — ; "^Rp is the azimuthal angle of the reaction plane. These angles are depicted in Fig. H] 
and in this figure the reaction plane is the plane of z = which is perpendicular to the 
direction of the magnetic field. 




FIG. 4: (color online). Schematic depiction of the transverse plane of a non-central heavy ion 
collision along the beam-axis (y-axis) (see also Ref. ^). The plane of z = is the reaction plane. 



In order to calculate the azimuthal charged-particle correlations, we follow the method 
of Ref. 121] to define the quantity A_|_ (A_) which is the positive (negative) electric charge 
difference in unit of e (— e) between on each side of the reaction plane, i.e., the z = plane 
in our notations. 

In Fig. [5] we give an schematic illustration of the chiral magnetic effect, and detailed 
discussions are presented in the caption. From this figure, one can clearly find that it the 
P-violating effect, i.e., the nonvanishing chiral chemical potential /is, that results in the 
difference of the numbers of the right-handed quarks (anti-quarks) and left-handed quarks 
(anti-quarks). Then, under an external magnetic field the number of quarks moving along 
the direction of the magnetic field (i.e., the number of the right-hand quarks in Fig. [5]) is 
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a (^5 = 0) 



FIG. 5: (color online). Schematic illustrations of the electric charge separation and the chiral 
magnetic effect (see also Ref. [12]). Here we take "up" quark u and its anti-quark u for example. 
The thick black arrows denote the directions of quark spins and the thin red arrows denote those 
of the momentum of quarks, a. Under the external magnetic field, quarks in the lowest Landau 
level (not including quarks in high order Landau levels) are polarized. The spins of u quarks are 
parallel to the direction of the magnetic field and those of u anti-parallel to that direction. In the 
case of /is = 0, i.e. these is no "P-violating effect, the number of right-handed quarks is equal to 
that of left-handed quarks and so the number of quarks above the z = plane is also equal to that 
of quarks below the z = plane. Therefore, when /X5 = there is no electric charge separation and 
chiral magnetic effect, b. In the case of /X5 7^ 0, the number of right-handed quarks is unequal to 
that of left-handed quarks and the numbers of quarks on two sides of the z = plane are different, 
resulting in the difference of the electric charges between on each side of the reaction plane, which 
is the electric charge separation effect. 

different from that of quarks moving against it (i.e., the number of the left-hand quarks in 
Fig- ED, and in this way the phenomenon of the electric charge separation takes place as the 
Fig. [5] b shows. 

Considering the simple system composed of only one type of fermion (with positive charge 
qe) and anti-fermion once more. From the Fig. O and the discussions above, we can easily 
find that the difference of the numbers of positive fermions on each side of the reaction 
plane is just the difference of the numbers of the right-handed fermions and the left-handed 
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fermions in the lowest Landau level, which is just 



5|n=0 



n=0 



71=0 



where 



R 



i + r 



ip and ipi 



1 - 7' 



(20) 



(21) 



2 ' ^"2 
We should emphasize that the subscript n = in Eq. (I2U]) indicates that the difference of 
the fermion numbers on the two sides of the reaction plane only comes from fermions in the 
lowest Landau level, since only fermions in the lowest Landau level are polarized, which is 
proved in Appendix |A] and Appendix [Bl The superscript + in Eq. ( 120|) means that only 
the positive fermions (not the negative anti-fermions) are included. Therefore, employing 
Eq. (1A46P in Appendix R] we can further express A^sI^^q 



5\n=0 



n=0 



27T ^ J 27r\ E ^ ' ^' 

LI ^ 



V 



\q\eB 



s=±i ■ 
00 



n=0 



P3 



P3 



1 



(22) 



where E is given by Eq. (IT7|) . In the same way, we can obtain the difference of the numbers 
of negative anti-fermions on each side of the reaction plane, i.e.. 



n=0 



q\eB 



-V 



2-n 

\q\eB 
47r2 



E 

s=±l 
00 



n=0 

rfp3 /s|e| 



27r 

P3 



E 



n=0 



'^^'Ee(^+/^-t^5)/T^i 



dp: 



P3 



.(23) 



Until now, we have obtained the difference of numbers of the positive fermions (negative 
anti-fermions) between on each side of the reaction plane, so the electric charge difference 
can be easily obtained as A+ = IgjA'sl^^Q and A_ = IgjA'sl'^Q. We should emphasize 
that though the electric charge differences A+ and A_ are the differences of quark electric 
charges in our picture, these electric charge differences are conserved through the hadroniza- 
tion processes (or other processes) and are observed in the heavy ion collision experiments. 
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because the hadronization processes (or other processes) are difficult to result in electric 

nn 

charge separations (for more discussions see Ref. |9|,ll0l). 

The calculations above can be easily extended to the 2+1 flavor quark system, and for 
this system we have 

eB r 



A, 



VNr 



+qs 



dp- 



dp- 



P3_ 

Ed 



i^d J^d 



P3 .M 



P3 



(24) 



and 



-VN, 



+c 



eB 
'4^ 



dp3^ f{E, - /X, - g-f,l) - f{E, - /z. + ^fif) 
dp3^ f{Es -lis- -p-Z^s) - f{Es- fis + l^/is) 



dp: 



P3_ 

Eu 



^3 .,d\ 



P3 



E 



(25) 



where Ef is given by Eq. ( HTl) with n = 0, and the distribution functions f{x) and f{x) for 
quarks and anti-quarks respectively, are given by Eqs. (HM ( |T5l) . One could find that when 
the quark chemical potentials are vanishing, i.e., /ij = (z = u,d,s), we have A+ = — A^. 
In the high temperature limit, the masses of quarks can be neglected and when the quark 
chemical potentials are vanishing, Eq. flM|) and Eq. fl25|) can be calculated analytically. The 
results are 



A_ 



-A_ 



f=u,d,s 



47r2 



V 



27r2 



(26) 



In Fig. [6] we show A+/V^, where V is the volume of the system, as function of the 
temperature in the PNJL model with eB = lO^MeV^, /Xj = 0, and several values of the 
chiral chemical potential /ig. In fact, the dependence of A^/V on the temperature is similar 
with that of the chiral electric current density. When the temperature is above the critical 
temperature of the chiral phase transition, approaches its high temperature limit 

value given in Eq. fl26|) . which is shown in Fig. [6] by the horizontal lines. However, once the 
temperature is decreased to that below the critical temperature, chiral symmetry is broken 
and quarks get large constituent masses, which results in that the electric charge difference 
between on each side of the reaction plane is suppressed drastically. 
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FIG. 6: (color online). Ratio of the positive electric charge difference between on each side of the 
reaction plane and the volume of the system, i.e., A-^/V, as function of the temperature in the 
PNJL model with eB = lO'^MeV^, /ij = (z = u,d,s), and several values of the chiral chemical 
potential yUs. The three horizontal lines denote the values of A^- /V in the high temperature massless 
limit, corresponding to fi^ = 150, 200, and 250 MeVs from bottom to top, respectively. 

In order to determine the azimuthal charged-particle correlations in heavy ion collisions, 
we need to calculate the quantity (A^-) which is the total positive (negative) electric 
charge number in unit of e (— e) on both sides of the reaction plane. Considering the simple 
system composed of only one type of fermion (with positive charge qe) and anti-fermion, we 
can easily find that the total positive (negative) electric charge is the sum of the positive 
(negative) electric charge of the right-handed and left-handed fermions (anti-fermions), i.e. 



q\i^J d^xipR-f^ipn + j d^xi>Ll^^L 
q\{^j d^xip'y'^ip^ I 

dp3 



2vr ^ ^ / 27r 

n=0 s=±l ' 



v\qf-f: 

27r ^ 



n=0 



dp3 



+ 



dps 



where we have 



2n\q\eB + pi 



(27) 



(28) 
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and 



E 



2n\q\eB + p| + m? 



(29) 



The superscript + on the right hand of the vertical hne in Eq. ( !27l) denotes that only 
fermions with positive charge are included (not including negative anti-fermions). We should 
emphasize that all Landau levels are summed in Eq. f l271) . which is different from the electric 
charge difference between on each side of the reaction plane in Eq. fl22l) . where only particles 
in the lowest Landau level contribute to the charge asymmetry. In the same way, one can 
also obtain 



N_ 



\Q\ 



n=0 s=±l 
oo 



dp3 / 
2n V 



2eB I" f dps 



v\q\'—y 



1 



277 p(£;+M-J#/^5)/T 



dps 



1 



Similarly, for the 2+1 flavor quark system we can obtain 



277 JE+ti+i^f,5)/T 



+ 1 



(30) 



n=0 s=±l 

, 2 f dp3 rfr. , s\ed\ a\ 



s\e 



Eu 



f(Es + Us 



(31) 



and 



n=0 s=±l ^ 



E. 



f {Es - fis- 



dp3 



(32) 



277 " ^ ^ E, 

We should comment that in Eqs. (!3T!) (!32|) we have assumed that the total positive (negative) 
electric charges of the quarks and anti-quarks in the fireball at early stage are conserved 
through the subsequent evolution of the QGP and are observed by the multiplicities of the 
produced charged particles in experiments. Although this is an assumption, it is reasonable. 
Because if the centrality is fixed and the collision energy is increased, on the one hand, the 
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temperature of the QGP at early stage is increased, which results in that the total positive or 
negative electric charges of quarks and anti-quarks increase, on the other hand, the increase 
of the collision energy will lead to the increase of the multiplicities of the produced charged 
particles. Therefore, the total positive (negative) electric charges of the charged particles 
produced in heavy ion collisions increase with those of the quarks and anti-quarks. 

So far, we can calculate the azimuthal charged particle correlations (cos(0q,-|-0^ — 2\l//jp)) 
in heavy ion collisions. With the notation aap = — (cos(0q, -|- 0/3 — 2\E'ijp)), it can be shown 
that 



i2i 

16 Ni ' 16 N': 



_ " \^+/ _ " \^-/ (oo\ 



and 



where the azimuthal angle distribution of the charged particles is assumed to be 

^ = T^iV± + 7 A± sin(0 - ^i^p) . (35) 

Since we mainly focus on the influence of the QCD phase transitions, especially the chiral 
phase transition, on the chiral magnetic effect in this work, we will neglect the screening 
suppression effect due to the final state interactions [l2] and make fii = {i = u, d, s), then 

we have a+_|_ = a = — a_|__. Therefore, we only study 0++ in the following. 

In Fig. [7] we show a_|_+ defined in Eq. ( l33l) as function of the temperature at several 
values of the chiral chemical potential /is (/is = /^^ = /x^ = /i^) and the magnetic field 
strength. We find that there is a pronounced cusp in a++ at the critical temperature during 
the chiral phase transition (the critical temperature = 209 MeV for /is = 150 MeV and 
Tc = 185 MeV for /is = 250 MeV in the PNJL model). From the Fig. [7] one can also find 
that although the value of a++ is proportional to the square of the magnetic field strength, 
the shape of the curve for function of temperature is almost independent of the 

magnetic field strength. Furthermore, the cusp at the critical temperature in the curve 
becomes much sharper with the increase of the chiral chemical potential. With the decrease 
of the temperature, when the temperature is below T^, chiral symmetry is dynamically 
broken and quarks obtain large constituent masses. We should emphasize that it is the large 
quark mass that results in the drastic suppression of the chiral electric current density (see 
Sec. lIIII) . the electric charge difference between on each side of the reaction plane (see Fig.|6]), 
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FIG. 7: Correlation as function of the temperature calculated in the PNJL model with fi^ = 
150 MeV (left panel) and /xs = 250 MeV (right panel). The magnetic field corresponds to eB = 
5 X 10^, 10^, and 5 x 10^ MeV^ from top to bottom, respectively. 

and the azimuthal charged particle correlations. Furthermore, the chiral magnetic effect is 



close related with the axial anomaly 12|. Without axial anomaly there would not be the 



chiral magnetic effect. Since the axial anomaly can be suppressed by the mass effect, which 



has been discussed in detail in Ref. 29|], it is natural to expect that the chiral magnetic 
effect can also be suppressed by large constituent quark masses. Therefore, it is reasonable 

that the azimuthal charged particle correlations described by a++ (a and a+_) defined in 

Eqs. ( j33l) (!34l) are quite decreased once the temperature is below the critical temperature. 
It can been seen from Fig. [7]that, when the temperature is above Tc, a++ decreases with the 
increase of the temperature, which is because higher temperature makes it more difficult to 
polarize quarks with magnetic field and thus suppresses the charge separation effect. 

What do our calculated results imply in future energy scanning experiments of heavy ion 
collisions? With the decrease of the heavy ion collision energy, the temperature of the QGP 
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produced in the fireball at early stage is also decreased. Since the magnetic field produced in 
„on-ce,...I con.io. decay, wit. .i.e ^ . ob..ved ch.ge .paction ,„.n„ ca.ne. the 
information of the QGP at early stage. Therefore, we expect that when the collision energy 
is decreased to a value that cannot drive the chiral phase transition, the azimuthal charged 
particle correlations (especially for the same charge correlations, because the opposite charge 
correlations are suppressed by final state interactions) are quite suppressed. So this property 
can be employed to search for where the QCD phase transitions take place. 

From the calculations of the azimuthal charged particle correlations (cos(0q, + 0/3 — 2\E'/jp)) 
above, we can find that there are some uncertainties on the total positive (negative) electric 
charges of the produced charged particles (A^-). This little defect motivate us to search 
for other better correlators which are not divided by the square of A^-|. or A^_. In fact, this 
kind of correlators has been proposed by D.E.Kharzeev and his collaborators I2I. In the 
following we compare these two kinds of correlators briefiy. First of all, we make (0^ ~ 
^/?p) 0a and (0/3 — "^Rp) — i- 0/3, and then the 0q, and 0^ are the azimuthal angles of 
produced particles with respective to the reaction plane as Fig. H] shows. For each collision 



event, we follow the definition of the correlators in Ref. 12], i.e.. 



^ No, 

fi^a, 0/3) = ^T^^ Yl Yl + (36) 

i=l j=l 

In the same way, here a, /3 = ± denotes the electric charge. In order to remove the multi- 
plicity fiuctuations the correlators are averaged over A'^g similar events. Then one obtain 

aa/3 = -(cos (0a + 0/3)) 
^ n=l 

where the correlators are those calculated in Fig. [71 Furthermore, correlators which 
are not divided by square of the total multiplicity of charged particles are also proposed by 
D.E.Kharzeev and his collaborators and are thought to be very useful, which are 



1 \ 

ba^ = -^Yg{(pa,(Pt3), (38) 



n=l 

with 



g{.<Pa, 'Pp) = COS(0„, + 0/3j) . (39) 

i=l j=l 
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Same as the aap, it can be found that 



b 



'+ 



■+ 



16 



16 



(40) 



and 



71' 



.2 



b 



16 



(A+A_). 



(41) 



In Fig. [6] we have calculated the A+/y as function of the temperature during the QCD 
phase transitions, and find that A+/y is rapidly suppressed and approaches zero when the 
temperature is below the chiral critical temperature. Within similar collision events (similar 
centrality and atomic number), it can be expected that the dependence of the volume of 
the fireball at the early stage on the collision energy is mild. Therefore, it can be predicted 
that with the decrease of the collision energy, the correlators ba/3, which are the azimuthal 
charged particle correlations not divided by the square of the total multiplicity of charged 
particles, get a sudden suppression at the critical temperature of the QCD phase transitions 
and approaches zero rapidly. 

V. SUMMARY AND DISCUSSIONS 

In this work, we have studied the influence of the QCD phase transitions on the chiral 
magnetic effect. The chiral electric current density, the electric charge difference between 
on each side of the reaction plane, and the azimuthal charged particle correlations in heavy 
ion collisions are calculated in the PNJL model, and their dependence on the temperature 
are studied. We find that with the decrease of the temperature, the chiral electric current 
density and the electric charge difference between on each side of the reaction plane are 
suppressed abruptly at the critical temperature of the QCD phase transitions and approach 
zero rapidly, since below the critical temperature the chiral symmetry is broken and quarks 
obtain large constituent mass. It is the large quark mass that suppresses not only the axial 
anomaly but also the chiral magnetic effect. For the azimuthal charged particle correlations, 
we study not only the correlators aai3, which are the correlators divided by the square of 
the total multiplicity of charged particles and are measured in current experiments, but also 
another kind of correlators bai3 which are not divided by the square of the total multiplicity. 
We find that both 0^/3 and 6^/3 get a sudden suppression at the critical temperature of 
the QCD phase transitions. Furthermore, the correlators 6q,^ approaches zero rapidly once 
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the temperature decreases to values that are below the critical temperature. Therefore, 
It indicates that azimuthal charged particle correlations (both a^^ and bap, in fact ba^ is 
better because the correlators ba/3 remove the fluctuations of the total multiplicity of charged 
particles) can be used as a signal to identify the chiral phase transition in the energy scan 
experiment in RHIC. 

We should discuss the possibility that using the azimuthal charged particle correlations 
(oq/j and bajs) to search for where the QCD phase transitions take place in future energy 
scanning experiments in RHIC. In order to simplify the calculations, we make the magnetic 
field strength and the chiral chemical potential fixed across the QCD phase transitions in this 
work. In the realistic situations the magnetic field decays with time and the chiral chemical 
potential has some distribution. However, we think that our simplification is reasonable and 
would not change our conclusions. The reasons are listed below: 

(1) In this work we are concentrated on the influence of the QCD phase transitions, es- 

pecially the chiral phase transition, on the chiral magnetic effect embodied by the 
phenomenon of the charge separation. Since the chiral phase transition takes place 
during a very narrow region of the temperature (or the collision energy equivalently) 
as Fig. [2] shows, the dependence of the magnetic field and the chiral chemical potential 
on the temperature is limited in this narrow region. 

(2) Indeed the magnetic field decays with time in a collision event, which has been confirmed 



in Ref. [1^. However, what influences on our calculations is the dependence of the 
magnetic field in the fireball at early stage on the collision energy (different events with 
different collision energy). Why is the magnetic field at the early stage of the evolution 
of the fireball? This is because the electric charge difference between on each side of the 
reaction plane A±. (describing the magnitude of the charge separation) is proportional 
to the magnetic field strength as Eqs. fl2^ fl25|) show (we should emphasize that the 
total multiplicity of the charged particle N± is almost not affected by the magnetic field 
with value eB = 10^ ~ 10^ MeV^ in the non-central heavy ion collisions). Therefore, 
the electric charge difference coming from the early stage of the fireball evolution is 
much larger than that from the late stage, since the magnetic field at early stage is 
larger than that at late stage. So we are more concerned about the magnetic field 
at early stage. Fortunately, for similar collision events (similar centrality and atomic 
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number), the dependence of the magnetic field on the colhsion energy is very mild 
(comparing Fig.A.l. (center of mass energy per nucleon pair being ^/s = 62 GeV) 
with Fig. A. 2 ( ^/s = 200 GeV) in Ref. |l2|). This is because what determines the 
magnetic field is the velocity of the heavy ion for similar collision events. However, 
for ^/s = 200 GeV the velocity of the heavy ion is f = 0.99995c where c is the light 
speed; for ^/s = 62 GeV the velocity of the heavy ion is f = 0.99948c. Therefore, 
although the difference of the two collision energy is quite large, the difference of 
their corresponding velocity of the heavy ion is quite small, which results in that the 
difference of the magnetic field is small. 

(3) As for the chiral chemical potential, we should comment that in the chiral symmetry 
broken phase, quarks get constituent mass. It is found that the mass always causes 
the asymmetry between the number of right-handed and left-handed fermions to de- 



cay 



31| . i.e., the chiral chemical potential decreases with time in the chiral symmetry 
broken phase. Our calculations above indicate that when the temperature is decreased 
and crosses the critical temperature of the chiral phase transition, large quark con- 
stituent mass suddenly suppresses the electric charge difference between on each side 
of the reaction plane. These calculations are performed with the chiral chemical po- 
tential fixed. If we further consider that the chiral chemical potential is reduced when 
the temperature is below the critical temperature, the suppression is much more sig- 
nificant. 

(4) In this work we perform our calculations with multi-values for the magnetic field 
strength and the chiral chemical potential (see Fig. |3l Fig. [6l and Fig. [7]). For all 
these values we find the same conclusion that in the chiral symmetry broken phase, 
the chiral magnetic effect is quite suppressed and almost vanishes. Therefore, the chiral 
magnetic effect can be used as an order parameter of the QCD phase transitions. 

(5) We should emphasize that the physical essence underlying our calculated results is very 

important. It is the large mass that suppresses the axial anomaly and the asymmetry 
between the number of right-handed and left-handed fermions, which is verified in 



general quantum field theory 29|, l3l|- When the mass of the fermion approaches 
infinity, there would be of course no difference between right- and left-handed fermions. 
Therefore, our calculated results is consistent with this basic principle and verify the 
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conjecture proposed by D.E.Kharzeev and his collaborators that the chiral magnetic 
effect can be used as an order parameter for the QCD phase transition 12 1. 
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Appendix A: Thermodynamics of a fermion system with V violation and under a 
background magnetic field 

Considering a system composed of only one type of fermion with positive charge e and 
mass m, and a homogeneous magnetic field with strength B is along the positive z direction. 
Assuming the system is in thermodynamical equilibrium with temperature T and chemical 
potential /i. In order to include the effects of V and CV violation, we follow the method of 



Ref. [15|] to introduce the chiral chemical potential /^s. We begin with the partition function 
of the system as 

z = T:^e-m-t^N-^.,N,)^ (Ai) 

where (3 = 1/T and quantities with hat are operators, and the Hamiltonian H is 

H = J d?xH = j d^xip{-ifdi - e^A' + m)ip, (A2) 

here i = 1,2, 3. The Hamiltonian density above can be obtained from the lagrangian 
density given by 

£ = ij{zj^D^ - m)ij, (A3) 

where = c?^ + ^eA^ and /i = 0, 1, 2, 3. In the Hamiltonian Eq. (]A2|) . we have used the 
fact that since we only consider the case with homogenous magnetic field along the positive 
z direction, we can assume Ai = —(1/2)5x2 and A2 = {l/2)Bxi. N and in the partition 
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function in Eq. (lAip are 



N 



(A4) 
(A5) 



respectively. 

First of all, we should solve the Dirac equation with a magnetic field, i.e. 



(A6) 



In the following, we employ the notations in Ref. 30|] and use the chiral representation of 
the 7 matrices, i.e. 



7 



1 

1 



Y 



a' 
-a' 



7 



-1 
1 



(A7) 



We express the four-component spinor as two two-component left-handed and right handed 
Weyl spinors, i.e. 

Then, for the positive energy solution, the Dirac equation (]A6P can be expressed as 



(A8) 



-id, - eA')a' 



m 



m 



E 



^1 



(A9) 



R 



We set = {—idi — eA'^)a^, then if we find an appropriate right handed spinor ipn, which 
is an eigenfunction of the with eigenvalue e, i.e. 



we have E"^ = e'^ + rn? and 



m 



(AlO) 



(All) 



We set ipR = y/E + e^"*, therefore, we have i/jl = VE — e^" and H^^'^ = e^"" 
For the negative energy solution of the Dirac equation flA6p . we have 



H_ m 

m —H_ 



tpR 



-E 



(A12) 



R 
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where we set = —H^ = {idi + eA^)a^ . In the same way, if we find a satisfying 
_ e^^('^)^ then we can obtain 



In the following, we will solve the eigenvalue equation 

if+r = er, (AM) 

and here, 

H+ = {-idi - eA')a' 

= -id^a^ + {-Ida - eA^)a'' 

= p^a^ + H^, (A15) 

where a = 1,2 and in the last line we have used the fact that is independent of X3 and 
A^ = 0, so the eigenstates in the X3 direction are free continuum of momentum. We should 
note that since {a^, H±} = 0, if there is a eigenstate |A) of H± with eigenvalue A > 0, there 
must be another eigenstate <j^\X) of H± with eigenvalue —A < 0. In the representation of 
I A) and cr^|A), Eq. ( ]A14I) can be expressed as 



where ^'^ = {ci,C2)^. Eq. flA16P has two solutions which are 



'"'^'^'^ 1. (A16) 



ci 1 ( sgn(p3)V|e| + A . 



and 



{-■ = ^1 -='5"<'")v^- I, with « = -N, (AI8) 



where |e| = a/A^ + p| and ^'^ has been normalized. Here we have use to stand for p^ 
without confusion. Equations (IA17P flAlSI) can also be unified to express as 



^ [s)sgn(p3)v/|e| + sA 
^ = , = , with e = s\e\. (A19) 

26 



In the same way, we can solve H_C^ ^'^^ = ^'^^ for the anti-particle, i.e., 



^-i, = ^^~^>^-(P^)^/^^^ I ,,th e = -.|e|. (A20) 



^2|e| V V|e| + sA 

Next, we turn to the eigen-equation of the transverse momentum H±\\) = A|A). It is 
obvious that we also have 

Hl\X) = X^\X). (A21) 

Hj^ can be directly calculated as 

The physical meanings of Eq. (1A22I) are very clear. The second line of Eq. ( 1A22I) indicates 
that the dynamics of particles in the transverse plane, which is perpendicular to the magnetic 
field, is the two dimensional homogeneous harmonic oscillation. The last line of Eq. flA22p 



includes contributions from the orbital and spin angular momentum in the z direction. 
Rescaling xi (ei?/2)^/^xi, X2 (ei?/2)^/^X2, andif| H\/eB, we find 

Hi = IM + X?) + (pI + xl)] - ik + a'), (A23) 

where Pa = —id/dxi {i = 1, 2) and I3 = XiP2 — X2Pi- In the following, we use the algebraic 
method to solve the problem of eigenvalue. Introducing annihilation and creation operators 

1 / 

tti = —={Xi + tpi), 

v2 

at = -^(x.-ip,). (A24) 

It can be easily verified that 

[ai,aj] = and [oj, a^] = (A25) 

thus these operators are annihilation and creation operators of boson. Employing these 
operators we obtain 

Hl = atai + ata2 + l- {h + (T^), (A26) 
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with 

I3 = i{—aia2 + a^ai). (A27) 

To diagonahze the orbital angular momentum l^, we introduce two another pairs of annihi- 
lation and creation operators, and the annihilation operators are 

a+ = - ^^2), 

a_ = -^{ai +ia2). (A28) 
v2 

In the same way, it can be verified that they are also bosonic operators. Consequently, we 
can obtain 

^3 = a!j]a+ — ala_, (A29) 

and 

a^oi + a^a2 = a^a+ + ala_. (A30) 

Therefore, we finally have 

Hi = 2ata_ + 1 - = . (A31) 

V 2(ala_ + 1) / 

Here d^a_ is the boson number operators, and its eigenvalue can be denoted as n (n=0,l,2...). 
Recovering eB we obtain the eigenvalue of Eq. (IA2ip with = 2neB. We should emphasize 
that from Eq. flA3ip one can find that states of n > are degenerate with two different spins, 
while for n = state, there is only one spin. Therefore, this means that only particles in 
the lowest level are polarized by the external magnetic field. 

Until now, we have solved the Dirac equation under a background magnetic field. We 
summarize the results here. The wave function of the fermion is given as 



with s = ±1. (A32) 



When n > 0, 



when n = 0, 




with e = s|e|; (A33) 



with e=|p3|, (A34) 

(0, 1)^ for p3 < 
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r-1 _ 



(0. 



^ (1,0)^ for p3<0 
The wave function of the anti-fermion is 



, ,-.,1)^ for »3 > 
r' = \ \. , with 6 = -b3. (A35) 



When n > 0, 



when n = 0, 



v'=\ ^ , with s = ±l. (A36) 



e^^ = ^= , , with e = -s|e|; (A37) 

f(0,l)^ for p3>0 
r^i) = <^ ' ' , with e = -b3|, (A38) 

[ (1, 0)^ for p3 < 

[ (0, 1)^ for p3 < 



(A39) 



In the equations above, we have 



A = V2neS, (A40) 



\^+pI = JineB + p|, (A41) 



E = v/|e|2 + m2 = ^j2neB +pl + m?. (A42) 

It should be emphasized that s = +1 represents the state that the spin of a fermion or anti- 
fermion parallels to its momentum, i.e. the right-handed state, while s = —1 corresponds 
to the state that the spin anti-parallels to the momentum and thus is the left-handed state. 
Therefore, the helicity of a fermion is s and that of a anti-fermion is —s. 

Employing the standard canonical quantization procedure, we can express the Hamilto- 
nian H in Eq. flA2p as 

H = yff:i: fpEiat-at-KK-) 

n=0 s=±l ^ 

= E E / pEiat^< + K^K - 1), (A43) 

n=0 s=±l ^ 

where af and 6f correspond to the annihilation operators of the fermion and anti-fermion, 
respectively. V is the volume of the system. Furthermore, in the presence of an external 
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magnetic field, we have 

dpi dp2 eB 




oo 
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n=0 



In the same way, we have 



n=0 s=±l 



and 



n=0 s=±l 

Substituting Eqs. (IA43P flA45p flA46p into Eq. flAip . after a simple calculation we obtain 

n=0 s=±l ^ 

/i5)/T] } + In { 1 + exp [ - + /i - ^/i5)/T] . (A47) 



Appendix B: Another approach with modified Lagrangian 

In the Appendix |X1 we have derived the partition function of a fermion system with V 
violation and under a background magnetic field. We will discuss this subject in another 
approach in this appendix. 

Absorbing the term including A^s in Eq. f lAip into the Hamiltonian if. We can obtain the 
modified Lagrangian density 

£ = - m + /i57%')^- (Bl) 

Then the Dirac equation in Eq. flA6p is also modified as 

id^^ = ^(-idi - cAOtV + 7°^ - t'/^s] ^. (B2) 

This Dirac equation can also be solved through the same method used in the above appendix, 
and here we just give the results. The wave function of the fermion is given as 

u^={ v/^^^^^^^ ^ . = ±1. (B3) 

VE + (e-/X5)e 
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When n > 

= - 

when n = 0, 

f (1,0)^ for p3 > 
e={ ' , with e=b3|, (B5) 

1 (0, 1)^ for p3 < 



~1 _ 



^(1,0)^ for p3<0 
The wave function of the anti-fermion is 



When n > 0, 



, (0, 1)^ for p3 > 
r' = r , with e = -\p,\. (B6) 

ri.O)^ for »^ < 



1 ^ (-s)sgn(p3)^|e| - s\ 

when n = 



e^^ = ^= , , with e = -se; (B8) 



In the equations above, we have 



r-''^*"-^': With .^-|p3l. (B9) 

t (1, 0)^ for p3 < 

I (0, 1)^ for p3 < 



The partition function is given by 

n=0 s=±l ^ 



(BIO) 



A = ^2neB, (Bll) 
|e| = ^A2+pi = ^2ne5 + pi, (B12) 
E = V(|e| -s/i5)' + m2. (B13) 



+ ln|l + exp [- (E + /i)/r]}^. (B14) 

We should emphasize that the expression of particle energy E in the above equation is given 
by Eq. 
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In the same way, we can obtain 



n=0 s=±l 



s e 



E 



'{s\t\ - /is) 
E 

s\t\ - /is) 
E 



Making the ensemble average of we have 



m 

eB 



(■5|e| - /Us) 



n=0 s=±l • 

[s|e|-/is) 1 



E e(^-/^)/^ + 1 
(s|e| - /is)" 



A/'s can also been directly obtained through A^s = Tdln Z /dn^. 
Next, we calculate the chiral electric current density |l5| : 



J3 = y I d'^xtPj'^ij 



-( 



EE 

n=0 s=±l 
s\e\ - fX5^^Sp3 
E 'Me 



dp3 
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)(if)r&: + ( 



. ^|e|-/is w^^ 



E 



Therefore, the ensemble average of js is 



is = (js) 
_ e^E 
27r 



EE 



dp3 
2n 



'S\e\ - IJ,5\/SP3 



( 



n=0 s=±l ' 

■s|e| - /^sx /SPs^ 



E |e| ^e(^+'^)/'^+ 1 ^ ^ 
Same as A^s, js can also been obtained through 



E |e| ^e(^-'^)/^ + 1 
1 /s|e| - /iswsps 



E ^Mer 



eT/91nZ+ 91nZ- 



J3 



where 



V dp3 

n LI ^ 



n=0 s=±l ■ 
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and 

1"^- = "^ EE /^(f + ln{l + exp[-(£H-ri/r]}). (B20) 

n=0 s=±l J \ / 

are the fermion part and anti-fermion part of the partition function in Eq. ( ]B14I) . respectively. 
The minus related with the anti-fermion part in Eq. (IB 19 1) is due to our convention in 
Eq. (!AT2l) . 

From Eq. fIBlSI) we can find that only the states with n = contribute to the chiral 
electric current density js, since the integral variable ps ranges from — oo to cxd for n > 0, 
while from — oo to or to oo for n = 0. Therefor, can be simplified as 

f dp3i^s\e\- fi5..sp3. 1 



E; ap3 .S t - fi5. /^P3\ 

s=±l 



E ^Mer e(^-^)/^ + 1 



( S\e\- ^I5 ^(SP3. 1 ( S\e\- H5 ^(SP3. 



e 



In order to understand the physical meanings of the several terms in Eq. ( ]B21I) . we just 
extract the term related with the positive fermion and calculate 



s=±l I I 







oo 

dp3 





P3 - /i5 1 



a/(P3 - /^s)^ + exp[(A/(p3 - /^s)^ + m2 - /i)/T] + 1 

+ / dp3 , ^^-^^ ^ 

J-oo V(P3 - /i5)^ + exp[(A/(p3 - /Us)^ + m2 - ^)/r] + 1 
= 0. (B22) 

From the above equation we can find that, in the modified Lagrangian approach, it is 
P3 — not p3 that judges whether a particle is right-handed or left-handed, i.e., if the sign 
of P3—fJ^5 is same as that of the particle spin along the z-direction, the particle is right-hand; 
if opposite, then the particle is left-hand. However, the direction of motion of the particle 
is governed by p3 not p3 — fi^. Therefore, the properties of right-handed or left-handed of a 
particle cannot determine the direction of motion of the particle. In another word, the left- 
right asymmetry cannot be correctly related with charge separation effect in this approach. 
So, we conclude that the modified Lagrangian approach is inappropriate to describe the 
chiral magnetic effect. We should point out that the approach given in Appendix [A] does 
not have this problem. 
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Next, we continue to calculate the chiral electric current density js in Eq. (lB2ip . The 
only nonvanishing contribution to js come from the last term in Eq. flB21l) . i.e. 



J3 



dps 



dp. 



lim A— >oo 



(B23) 



27r2 

This is the result obtained in Ref. |15|. From our calculations above, one can find that ja 
calculated in the modified Lagrangian approach does not comes from the finite temperature 
contributions, however its source is the Dirac Sea which can be clearly seen in Eq. f lBlTp . 
Therefore, it is argued in Ref. 15] that is independent of the temperature, chemical 
potential /x, and even the mass of the particle (from the calculations above we find that 
the independence of js on the mass of the particle is due to the fact that the ultraviolet 
momentum in the Dirac Sea makes the mass of particle meaningless). We should comment 
that the results obtained in the modified Lagrangian approach are contrary to our physical 
intuition. On the one hand, the physical observable (here is the js) is unlikely to come from 
the Dirac sea. On the other hand, the chiral electric current density is also unlikely to be 
independent of the properties of current carriers. In this work we will show that the chiral 
electric density js calculated in an appropriate approach, i.e., the approach in Appendix \K\ 
comes from the finite temperature contributions not the Dirac Sea, and js is also dependent 
of the properties of current carriers and is also influenced by the external environment. 
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